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Abstract 

We examine here some of the effects on planetary 
spectra that would be produced by departures from isotropic 
scattering. The phase function w(l + a cosO) is the 
simplest departure to handle analytically and the only 
phase function, other than the isotropic one, that can be 
incorporated into a Chandrasekhar first approximation. This 
approach has the advantage of illustrating trends resulting 
from anisotropies while retaining the simplicity that yield? 
physical insight. 

An algebraic solution to the two sets of anisotropic H 
functions is developed in the appendix. It is readily 
adaptable to programmable desk calculators and gives 
emergent intensities accurate to 0.3 percent, which is 
sufficient even for spectroscopic analysis . 
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I . INTRODUCTION 

A considerable literature exists on the photometric and 
polarization effects of anisotropic scattering by planets. 
But spectroscopic effects have been examined quantitatively 
almost entirely with isotropic-scattering models. An 
important exception is the work of van de Hulst and 
Grossman (1968), who initiated a large number of numerical 
solutions with the doubling method and the Henyey-Greenstein 
phase function. 

The problem is not that spectroscopic data are more 
difficult to analyze with anisotropic scattering. Rather, 
the "trouble" is that, virtually any reasonable model of 
atmospheric inhomogeneity and scattering phase function can 
be made to duplicate the data well enough, although one's 
conclusions may be dependent on the adopted model. The 
trade-off of one parameter for another allowed van de Hulst 
and Grossman to develop their "similarity relations". But 
they noted, "As a consequence, it will be difficult, in 
practice, to interpret absorption lines observed in diffuse 
reflection in a unique manner." 

The investigation of other analytic models may help to 
illustrate weaknesses in interpretations that are made on 
the basis of isotropic, homogeneous theory. The phase 
function £5(1 + a cos 0) is not very representative of Mie 
scattering with a strong forward lobe, such as we find in 
Venus’ atmosphere. Nevertheless, it will illustrate some of 


the effects to be expected (e.g., the role that anisotropy 
of scattering can have in causing absorption bands to vary 
with position on the disk and the planetary phase) . 


II. REFLECTED INTENSITY FOR THE ANISOTROPIC 
PHASE FUNCTION: A SUMMARY 


The radiative transfer problem with a phase function 
for single scattering of 65(1 + a cos©) was solved exactly 
more than 25 years ago (Chandrasekhar, 1950 - hereinafter 
referred to as "RT") . The solution for the emergent intensity, 
I, is divided into two parts: one, I^, independent of the 

azimuth of the incident sunlight (relative to the observer) , 
and the other, I dependent on azimuth. Each of the 
components is solvable in terms of particular Chandrasekhar 
H- functions. Thus, having tabulations of these H-functions 
is tantamount to having numerical solutions for any given 
parameters, 65 and a. 

In the conventional notation (RT) , the equation of 
transfer for I(t;u,iJj) is 


^ dT 4 IT 


J 


'1 C2ir 


■1 -J 0 


I p(y,y'; tjj-ifj' ) dfdy ' 


( 1 ) 


_ F -x/y 

T a 


o P(Ur-V 0 ; ^-^ 0 ) 



With a scattering phase function. 


p(y,y'? = 5(1 + a cos0) , -1 < a < 1, (2) 

and 

cos0 = y y' + (1-y 2 ) (1-y ' 2 ) '^‘cos (ip-il* 1 ) , (3) 

we write 

I (t;y ,tf>) = I*°*(T?y) + I^(T;y) cos (ij/-^) . (4) 


Equation (X) now yields separate equations for the 
azimuth-independent intensity, I^, and the azimuth-dependent 
amplitude, 1^ . Both equations may be solved in the first 
approximation in a manner analogous to, but rather more messy 
than, the solution for isotropic scattering (RT, pp 149ff) . 

The exact solutions for the emergent intensity have 
also been obtained. For I^O) W e have 


I (0) (0 , y ) - H (0) (y)H (0) (y Q ) [l-c o (y+y o )~a(l-w)yM o ], 


(5) 


where (RT, p 141) 


c o 


5 a j 0 ^ a ( 1-65 ) 


2-w a 


( 0 ) 


( 6 ) 


( 0 ) 


and where a vu/ is the nth moment of (y). 


n 
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The solution for 1^ is 


l (1> (0,y) - JJiL. {x- W *> 1/2 U-M/) ' /2 H (1) <„_) . 


(7) 


The integral equation for H ^ (y ) (where j = 0 or 1) is 


H (y) 


f 1 

, 1-2 
■i — -J 0 


H' (j) (y • )dy 


>/■ 


y’V (j) (y ’ )H (j) (y') 


y + y' 


dy*. (8) 


(j) 


Here the "characteristic function", r J/ , is related to kj , 
the root of the "characteristic equation" that gives the T 
dependence of the equation of transfer (1) , by (RT, p 114) 


(j) 


1 . 2 , , 2 \ 1 / 2 


V SJ ' (y) = i a-k. z y 2 ) 


We have in the first approximation 


k o - 57 t(l-ffia W /)] ,/2 . 


and 


k n = — 

1 y. 


wa(l-y 2 ) 
1 __J — 


1/2 


where y ^ = 1//3 is the first Gaussian division point. Hence 
if we write the characteristic function as 


(y ^ ) = rj + s j U j 2 ' 


( 9 ) 


( 10 ) 


( 11 ) 


we obtain r Q = 2/2, s q = ua(l-w)/2, and -3^= 5a/4 


( 12 ) 
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In Tables 1 and 2 we give solutions for H ^ and its 
moments obtained by successive evaluations of equation (8). 

The solution was iterated until successive values of the 
zeroth moment, a^, converged to 1 part in 10 4 with a 
seven-point Gaussian quadrature. When used with Chandrasekhar's 
tabulations of for a = 0 (RT, p 125) , a = +1 (RT, p 139) , 
their moments (RT, pp 126, 328, 141), and (RT, p 141), 

these tables provide a convenient basis for tabulating solu- 
tions to equation (4). Additional values of H ^ f or hi = 0.975 
and several values of a between -1 and +1 were tabulated by 
Harris (1957) . 


III. PLANETARY SPECTRA 

In addition to evaluating solutions to the integral 
equations, I have developed (cf. Appendix) an approximate 
analytic solution to equation (4), which is accurate to 0.3 
percent at worst. With this expression I have made a number 
of sample comparisons of planetary spectra with strongly 
forward (a = +1) and backward (a = -1) scattering functions. 

Let 0 be the polar angle (co-latitude) of the observed 
point, in a coordinate system whose pole is the pole of the 
ecliptic. Also let <J> be the longitude of the point measured 
from the sub-earth point; 4> is the planetary phase angle. 

Then x and X Q are the zenith angles (at the point of observa- 
tion) of the earth and sun. That is. 


i 

( 


p s cos X - sin 6 cos <J> , 


(13) 


P = cos x o = sin 6 cos (<!>-<{>) , 


(14) 


and the azimuthal angle of scattering is given by 


cos (t|/-4> 0 ) 


PPo - cos4> 

-lyr 


iJT 


(1-p ) ' (1“M„ ) 


(15) 


With this geometry we obtain the spectral absorption 
curves in Fig. 1 for a weak line. The corresponding curves 
for isotropic scattering (a = 0) are not drawn in, but they 
lie almost mid-way between the extreme anisotropic cases. 

Generally speaking, the limb/center ratios are not 
strongly dependent on the latitudinal zone but are sensitive 
to the anisotropy factor. At 60° the limb/center absorption 
ratio is 0.50 for backward scattering (a = -1) and 0.70 for 
a = +1. At the crescent, there is, of course, much less 
variation, but the trend is reversed: ~ 1.0 for a = -1 and 

0.92 for a = +1. 

Similarly, we can make some useful estimates of how the 
phase variation will change with a. For weak lines a simple 
approximation, treating the continuum reflection as weighted 
over the planet like a lambertian reflector, predicts that 
every zonal band has the same phase variation (Chamberlain 
19 75) . This uniformity, in fact, is nearly the case in our 
model calculations, but again the variation is strongly 



dependent on a. The phase absorption ratio, computed for 
(<£ = 60 0 )/C<t> = 150°) is 1.5 for a = -1; 2.25 for a = 0; and 
3.3 for a = +1. 

That the spectral phase curve will go through a 
secondary minimum near full phase seems unlikely, at least 
with the simple anisotropic phase function examined here. 

For purposes of comparison with the complete phase curves 
for isotropic scattering (Chamberlain 1970, Fig 6), the 
albedos adopted in the present example (w„ = 0.975, = 

0.99) correspond to the parameters used in that paper of 
w = 0.0154 (very weak line) and q = 1.54 (moderately important 
continuum absorption) . 

Finally, the curve of growth for continuum albedos ft = 1, 

~ 1 / 2 

is very nearly linear with (1-tu^) , for observations at 

any one point, as it is in the isotropic case (Chamberlain, 
1965) . However, the slope varies with the anisotropy factor 
and to some degree with the scattering azimuth. 

Further discussion of the calculations will be more 
illuminating in the context of reliable data or comparisons 

with more realistic (and complicated) phase functions. It 
will also be necessary to consider what effects vertical 
inhomogeneity will have on the spectral phase variation. It 
is not clear that the effects of inhomogeneity and anisotropy 
can be distinguished. We expect to report on these matters 
at a later date. 
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APPENDIX: ALGEBRAIC SOLUTIONS 

?OR H-FUNCTIONS 

In this Appendix we develop practical formulae for 
computing the H-f unctions and intensities quickly and 
accurately, without need for interpolation. The algebraic 
solution here is a generalization of the Chandrasekhar 
"first approximation", which can handle phase functions 
expressed as Legendre polynomials of order N=0 and 1. To go 
to order N requires (RT, p 151) an nth approximation given 
by n > N/2 + 1/4 . Being restricted by choice to the first 
approximation is equivalent to limiting ')'(y) to the two 
terms given in equation (12) . 

The technique developed here was previously introduced 
for isotropic scattering (N=0) in a paper concerned with 
coupled spectral lines (Chamberlain and Smith, 1975) . In 
that problem the parameters y and y Q are no longer simple 
direction cosines and may assume values exceeding unity. 

The present technique has the useful quality of providing 
H-f unctions quite accurately for large y and y Q , whereas 
linear interpolation from a table of H(y) can in such cases 
be very inaccurate (Chamberlain and Wallace, 1974, Appendix 
B) . 

The first approximation, as obtained from a one-point 
Gaussian summation, is 

1 + y/Vj 

nnqir ' 


H (j) (y) 


(A.l) 
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where j is an index to signify either the H ^ - or 
H^ - functions of equations (5) and (7). However, is not 
an especially accurate weighting factor over the whole range 
of y. Hence we could get a better representation by 


H (j) (y) 


1 + A^ (to, a) y 
1 + k..y 


(A. 2) 


and evaluate A to give the best fit to some specified set 
of points. 

However, with a little additional effort we may obtain 
a "second approximation" by putting equation (A. 2) under the 

integral of equation (8) . It is then advantageous to 

( i ) . . . 

evaluate A J after the integration, in order to optimize a 

fit with the second approximation. 

Although a bit messy, the integration is straight- 
forward . We have 


H (j) (pi 


l-2r . ^ 

D 3 


2s^| 1/2 , C 2 

+ ? o + JcT 


+ \ 


|ry 2 -.?y\ / S a /k.-S 2 /k. 2 

\-t=i 4r) ¥ + 1 A.p - 3 - j ' 


(A. 3) 


where 




s . s 


.A^ (w,a) s . A ^ 3 ^ (to , a ) 


3k 


_ _JL 


2k r* 
3 


{1 + k^y) 


(A. 4) 


s .y s . A ^ 3 ^ (S,a) y 




(1 + k^y) , 


(A. 5) 
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and 


S 


2 


Tj A^ (£3, a) 



s . A^ ) (£3, a) y 

■ J i. 


k. 

3 


1 - 


(1+k^y) 

k i» 


The r j and Sj are given below equation (12) . 

Optimum values of A^ (5, a) have been obtained by a 
trial-and-error adjustment of A^ and a fit of the resulting 
moments a^(ai,a) to the exact values for a=0 and +1 and « = 0, 
0.4, 0.8, and 1.0. Writing 

A (0) (w,a) = A 0 ll 4- B o a) , 

we obtain 


A - A (0) (Si, a) = /3 + 0.1853 £3 - 0.1251 w 2 4- 0.1851 £5 3 


and 


B, 


= - 1 = -0.0210 S + 0.0041 w 2 + 0.0169 £u 3 . 

A. (u>,0) 


The above three equations imply that A^ (w,a) is 
symmetrical in ± a and, indeed, to the accuracy of this 
approximation, that is the case. However, since a appears 
intrinsically in k , neither H^(1J) nor is symmetrical 


in a. 


To summarize, the algebraic solution for I 


(0) 


consists 


of A^(£3,a), computed from the above three equations, and 


H 


(0) 


(y) from equations (A. 3) through (A. 6) 


Here k is 
o 


given by equation (10) with y^= 1//3, and r o and s q are 


(A. 6) 


(A. 7) 


(A. 8) 


(A. 9) 
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defined by equation (12). Then c Q is given by equation (6) 
with and obtained by Gaussian quadratures of. 

H^(y). Finally 1^(0, U) is computed from equation (5). 
A useful check on the accuracy of the calculations is 
provided by the relation (RT, p.109) 


1 




+ 



) = 0 . 


We obtain A^ (wa) by fitting the zeroth moments a ^ to 

the moments of the (y) -functions computed from equation 

( 1 ) 

i8) . Values of H (y) have been tabulated by Chandrasekhar 
(RT, p 141). Additional values of H^(y) and its moments 
have been computed by the author. We find, over the range 
-1 < wa < 1, 

A (1) (5a) = /3 + 0.1055 (wa) + 0.0143 (wa) 2 . 

Thus l^(0,y) is found from equations ( 7) , (11) , (A. 3) to 
(A. 6) , and (A. 11) , 

In the worst cases has absolute accuracy to about 

3 units in the third decimal. This error is due to excursions 
of H^(y) about the precise values for a given w and a. The 
values of the moment a Q are good to 1 or 2 units in the 
fourth decimal. Thus the principal errors arise from the 
algebraic representation of equation (A. 3) and not from the 
formulae for A(w,a), equations (A.7)-(A.9). The H ^ functions 


(A. 10) 


(A. 11) 


are even slightly more accurate, so that the limiting accuracy 
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of I(y,y Q ) depends on . 

As noted at the beginning of the Appendix, in problems 
involving spectral lines that are coupled by having a common 
upper level, we concounter functions H(x) , where x is often 
greater than unity. The solution given here is especially 
useful for this situation. In fact, in the limit x •* 00 , we 
find 

H (x) *► (1 - 2r - | s)" 1 /*. {A. 12) 

This limit is obvious from equation (8) , but it is equally 
valid for equation (A. 3). [To prove the limit in the latter 
case it is necessary to expand both ln(l + 1/y) and (1 - 
1/ky)" 1 to order y" 2 .] The accuracy of the algebraic 
approximation for large arguments is very good, as shown for 
isotropic scattering by the Figure in the Appendix of 
Chamberlain and Smith (1975). 
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TABLE 1 


H (0) (u) FUNCTIONS COMMUTED FROM THE EXACT 
INTEGRAL EQUATION 


(a « -1.0) 


• 

- 

u 

, 0-2 

0.4 

wm 

0.8 

T 

0.9 

0.925 

0.95 

0.975 

o | 

1.0000 

1.0000 

i 

1.0000 

1.0000 

1.0000 

1.0000 

1 1.0000 

1.0000 

0.1 

1.0217 

1.0482 

1.0819 

1.1291 

1.1638 

1.1752 

1.1888 

1.2063 

.2 

1.0318 

;• .0726 

1.1273 

1.2092 

1 1.2742 

1.2965 

1.32 35 

1.3598 

.3 

1.0387 

1.0896 

1.1605 

1.2718 

1.3650 

1.3979 

1.4387 

1.4947 

.4 

1.0437 11.1025 

* 

1.1865 

1.3235 

1.4431 

1.4864 

1.5409 

1.6175 

.5 

1.0477 : 

>1.1128 ; 

1.2077 

; 1.3673 

1.5116 

1.5650 

1.6333 

1 . / 309 

.6 

1.0508 

1.1212 

1.2255 1 

1.4053 

'1.5727 

1.6359 

' 1.7176 

1.8365 

.7 

1.0535 

1.1282 

1.2406 

1.4385 

1.6276 

1.7002 

1.7952 

1.9:55 

.8 

1.0557 

1.1342 

1.2537 

1 . 4o80 

1 6773 

1.7591 

1.8670 

2.0287 

.9 

1.0575 

1.1394 

1.2651 

1.4943 

‘ 1.7227 

1.8133 

1.9338 

2.116'’ 

0.95, 

1.0584 

1.1417 

1.2703 ; 

1.5064 

1.7441 

1.8387 

1.9654 

• 

2.1589 

1.00 

I .0592 

1.1439 

1.2752 

1 1.5180 

' 1.7644 

1.8632 

1.9960 

2.2001 





(a * -0 

.5) 




— ^ — 

l.UUUT 

! l .oouO 

1.0000 

1 .OO 60 

| 1.0000 

1.0000 

| l.OOro 

1 . 00 O 0 

0.1 

1.0236 

j 1.0517 

1.0865 

1.1337 

1.1679 

1.1788 

1.1919 

1.2086 

.2 

1.0354 

1.0790 

1.1360 

; 1.2185 

1.2825 

1.3041 

1.3301 

1.3648 

.3 

1.0434 

1.0985 

1.1727 

1 

1.2856 

1.3776 

1.4097 

1.4492 

1.5029 

.4 

1.0495 

1.1136 

1.2020 

1.3414 

1 1.4600 

1.5023 

1.5553 

1.6289 
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im 

U 

. 

0 . *4 

i 

0.6 

0.8 


*T* ' ’ ■— 

0.925 

0.95 

• . 

• ^ 

1.0543 

1.1256 

1.2260 

1.3892 

1.5328 

1.5852 

1 6517 

1.7458 

.6 

1.0583 

1.1356 

1.2*46*4 

1.2637 

1.4308 

1.5980 

1.6602 

1.7400 

1.8550 

• 

1.0615 

1.1*4*40 

1.4674 

1.6569 

1.7286 

1.8217 

1.9577 

• » 

1.0643 

1.1512 

1.2789 

1.5001 

1.7105 

1.7914 

1.8975 

2.0547 

.5 

1.0667 

1.1575 

1.2922 

1.5294 

1.7596 

1.8495 

1.9682 

2.1465 

.95 

1.0677 

1.1603 

' 1.2982 

1.5429 

1.7827 

1.8768 

2.0018 

2.1906 

1.00 

1.0688 

1.1630 


K5H 


-A.-2Q2L 

. 2.0344 

2-2.236 





(a « +0.5) 




0 


|H| 



1.0000 

1.0000 

1.0000 


0.1 

1.0276 

1.0592 

1.0965 

1.1442 

1.1768 

1.1870 

1.1989 

1.2138 

.2 

1 .0426 

1.0929 

1.1550 

1.2396 

1.3012 

1.3213 

1.3452 

1.3763 

.3 

1.0533 

1.1177 

1.2000 

1.3170 

1.4064 

1.4365 

1.4729 

1.5214 

.4 

1.0616 

1.137*4 

1.2366 

1.3827 

1.4990 

1.5391 

1.5883 

1.6551 

.5 

1.0682 

1.1536 

1.2673 

1.4398 

1.5819 

1.6320 

1.6942 

1.7802 

.6 

1.0737 

1.1671 

1.2936 

1.4903 

' 1.6571 

1.7170 

1.7922 

1.8979 

.7 

1 .078*4 

1.1787 

1.316*4 

1.5354 

1.7258 

1.7953 

1.8836 

2.0094 

.8 

1.082*4 

1.1887 

1.3365 

1.5759 

1.7890 

1.8680 

1.9692 

2.1154 

.9 

1.0859 

• 1.1975 

1.35*43 

1.6126 

1.8473 

1.9356 

2.0497 

2.2164 

0-95 

1.087*4 

1.2015 

1.3624 

1.6297 

1.8749 

1.9677 

2.0881 

2.2652 

L.O 

1.0889 

i 1.2052 

1.3702 

1.6461 

' 1-9015 

1.9987 

2.1255 

2.3128 






TABLE 2 

Moments of H^ U; (y) and the Constant 
(a - -1.0) 


u 

a<°> 

0 

»i 0> 

«< 0) 

c 

0 

0 .? 

l. 0434 

t 

0.5255 

0.3512 

-r 0469 

.4 

1.1027 

0.5690 

0.3764 

-0.0864 

.6 

1.18S7 

0.6140 

0.4144 

-0.1146 

.8 

1.3381 

0.7071 

0.4820 

-0.1217 

.9 

1.4758 

0.7962 

0.5476 

-0.1067 

0.925 

1.5279 

0.8305 

0.5732 

-0.0982 

0.95 

1.5955 

0.8755 

0.6068 

-0.0359 

0.975 

1.6943 

0.9421 

0.6571 

-0.0660 



(a * -0 

.5) 


0.2 

1.0495 

0.5293 

0.35^0 

-0.0237 

.4 

1.1145 

0.5683 

• 0.3818 

-0.0439 

.6 

1.2067 

0.6248 

0.4223 

-0.0588 

.8 

1.3591 

0.7208 

0.4921 

-0.0632 

.9 

1.4965 

0.8100 

0.5580 

-0.0558 

.925 

1.5479 

0.8439 

0.5833 

-0.0515 

• 95 

1.6140 

0.8880 

0.6163 

-0.0452 

0.975 

1.7099 

0.9528 

0.6652 

-0.0349 
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Ca = +0 

.5) 


A* 

U 

e 

O 

o 
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Fig. 1 - Variation of absorption across the visible disk 
of a planet at two phases. The "mid-latitude" curves apply to 
a zonal band at latitude 45°. The absorption s = 1-1 /I , is 

V V C 

for a weak line: specifically, = 0.95 with weakly absorbing 

tontinuum, = 0.99. 
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